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ABSTRACT 


This report covers the period of the contract extension, fron 
1 September 1979 to 31 August 1980. Tasks established and carried out 
during this period were: the development of a general ray trace evalua- 
tion program, called ICARUS; the study of novel Fresnel concepts; and 
the review of a report draft on novel Fresnel concepts. ICARUS is 
documented herein; reports on the novel Fresnel concepts have been 
submitted earlier. 


NEW TECHNOLOGY 

Fresnel Computation Program 

A general ray trace evaluation program, designated ICARUS, has 
been developed, with the goal of widening the generality over that in 
commercially available optical evaluation programs. The main thrust of 
the effort has been to expand the number of geometrical shapes, tilts, 
and decenters that could be explored, to include conic, aspheric, and 
toric surfaces. To keep the program task within time and budget, the 
effect of facet errors was not to be included. This task has been per- 
formed principally by Steve Eckhardt, graduate student, under the super- 
vision of A. B. Meinel. Documentation is appended to this report 
(Appendix A) . 

A review of the program led to the decision by us and by JPL's 
program management to terminate this ta'k before completion; however, 
we were able to bring it to a usable state; examples of Fresnels that 
can be ray traced are shown in Appendix B. ftie reason for withdrawing 


from the hope of a fully general program was the conclusion that a spe« 
cific program could be written for any novel Fresnel surface that could 
not be handled by an existing program. 

It is important to know what has not been included in the program 
to date. Not included are the following capabilities: (1) finite groove 
size* (2) diffraction, (3) local groove effects (facet tilt error), 

(4) substrate effects (lens surface errors) , and (5) general aspheric 
substrate shapes. The facet tilt error and surface errors are relatively 
easy to handle for reflective Fresnels in that whatever the angular dis- 
tribution of assumed errors is, the effect at the focus is twice the 
angular distribution, subject to the gecmietrical cosine effect of 
oblique incidence on the focal plane. The effect is considerably more 
ccnnplicated for refractive Fresnels. We will examine the analytical ex- 
pressions in future work if JPL requests it. 

Questions of adding facet error statistics have yet to be incor- 
porated. The analytical approaches to this are currently being evalu- 
ated and could appear as an addendum at the appropriate time. 

Fresnel Optic Study 

The rest of the cxirrent funding was expended in addressing novel 
Fresnel concepts, on which reports have been submitted, and in re- 
viewing a report draft on this topic prepared by JPL. Several of the 
novel concepts appear to be patentable, and a report was written on this 
question, but no reply was received concerning JPL's desire to proceed 
with this question. In the absence of a reply we will, at an appropri- 
ate time, prepare a paper for journal publication on these novel concepts. 
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APPENDIX A 


Prograa Title: 

Progran Code Naae: 
Writer: 

Date of Completion: 
Source Language: 
Availability: 

Abstract : 

Method of Solution: 
Program Capabilities: 


PROGRAM DOCUMENTATION 


PROGRAM IDENTIFICATION 


C^tical Sciences Center Fresnel Evaluation 
Program 

ICARUS 

Steve Eckhardt 
Optical Sciences Center 
University of Arizona 
Tucson, Arizona 85721 

May 1980 

Fortran (Data General FORTTIAN 5) 

Program listing as computer printout (attached) 

ICARUS is a ray trace program of greater gen- 
erality than ccmanercially available programs. 


PROGRAM DESCRIPTION 


Paraxial rays are traced by normal methods. 

Real rays are traced in accordance with the 
conic ray trace equations given in Table 1. 
Fresnels are treated by asstiming that the point 
of intersection of the ray with the surface is 
in the tangent plane of the lens and that the 
angle of the nozval to the ray is the same as 
if the Fresnel were an ordinary conic. 

Paraxial and real ray tracing, as well as 
radiant energy distribution computation, for 
conic, fresnel, toric, curved fresnel, toric 
fresnel, and curved toric fresnel surfaces. 

Not included are finite groove size, diffrac- 
tion, local groove effects (facet tilt error), 
substrate effects (lens surface errors) , and 
general aspheric substrate shapes. 
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Data Inputs: 

Program Options: 
Printed Output: 
Other Outputs: 

Flow Chart: 
Sample Run: 


Computer Equipment: 
Peripheral Equipment: 
Source Program: 

Variables and 
Subroutines: 

Data Structures: 

Storage Requirements: 

Maintenance and 
Updates: 


Data are entered interactively at the terminal, 
in response to quest icms posed by the program. 
Tables 3 and 4 provide guides to input options 
and format. 

See Table 3 

N/A 

A file is created for each new lens. These can 
only be deleted outside of the program. 

Figtire 1 is a general system flow chart. 

Table 3 is an annotated "walk-through" of the 
program. 


SYSTEM DOCUMENTATION 


Data General Eclipse 
N/A 

Table 5 is a complete source code listing for 
ICARUS and its subroutines. 

See the source code listing for variables. 
Table 2 is a list of subroutines, with descrip 
tions of their functions. 

N/A 

Minimal 

None 
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Fig. 1. ICARUS System Flow Chart 
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Fig. 1. ICARUS System Flow Chart (Continued). 
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Table 1. Ray Trace Equations for Conics. 


Opening ; 

A - [t2 ♦ (hy - . (hjj - 



N . t/A 

Transfer to Surface : 

B « c(x^^ + y^^) 

C = N - c(Lx^ + 

2 2 U 7 

D * [C - cB(l + kN )]^ D <0; ray missed surface 

A * B/(C + D) 

Xg » ♦ AL 

yg * y^ + AM 

Zj = AN 

Refraction : 

E - [1 - 2ck25]'** 

cos I > DE 

w * n/n’ 

2 2 2 W 

6 * (1 - p ♦ V cos I)^ - ycosi 

1 ■> ♦ y^cos^l < 0; 

total internal re f 1 ect i on 
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Table 1. Ray Trace Equations for Conics (Continued). 


a ■ -cE Xj 

P • -cE Xg 

Y • (1 - C(k + l)Zs)E 

L' ■ vL ♦ 6o 

M* > vM 50 

N' « mN ♦ 6y 


Transfer to Tangent Plane ; 

Xt * Xg ♦ L(t - ig)/N 
y^ - Xg + M(t - Zg)/N 

A « [(t-Zs)2 ♦ ♦ CXfXg)2]*s 

Variable Definitions ; 

ray height in x on tangent plane 
X( ^ ray height in y on tangent plane 

Xg « ray height in x on surface plane 

Xg « ray height in y on surface plane 

t ■ distance from object to first lens surface, or distance 

between surfaces 

\ ■ object height in x 

hy ■ object height in y 

L,M,N = x,y,z direction cosines for the ray before refraction 

L*,M*,N' « x,y,z direction cosines for the ray after refraction 

a,0,Y ■ direction cosines of the normal to the surface in x, y, 
and z, respectively 

A B ge(»etrical path length between a surface and the following 
tangent plane; also vice versa, the refractive index times 
the optical path length 

I ■ the angle of incidence on the surface 

n,n* B the refractive indices of the materials before and after 
the surface 

B,C,D,E,5 > dummy variables 
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Table 2. 

List of Sid}routines. 

Name 

Called by 

Description 

CHANGE 

Main progran 

Interacts with user to make individual 
changes in the lens without requiring re* 
typing of entire lens; inputs are as in 
subroutine SltXlE 

LENR 

(lens read) 

Main prograa, 
PARAX 

Reads a currently existing lens file into 
the active area 

LENT 

(lens type) 

Main prograa 

Displays the current lens on the terminal 
in tabular form 

LENW 

(lens write) 

CHANGE, 

STORE 

Creates a file called by any name the user 
specifies, and writes the lens to it. 
Should the need arise, this file nay be 
accessed for modification without use of 
the program 

PARAX 

Main prograa 

Does a paraxial ray trace that enables the 
user to locate the Gaussian image plane, 
and in general to do a rough lens design 

RAY 

Main prograa, 
RED 

Traces a specified ray through the lens 
(see TORT and SYI40 

RED 

Main program 

Calculates the radial energy distribution 
in the image plane 

SORT 

RED 

Orders the rays by increasing radial dis- 
tance from the paraxial image point 

STORE 

Main prograa 

Interactively accepts data input 

SYKM 

RAY 

Does refraction calculations for symmetri- 
cal surfaces 

TORT 

RAY 

Does refraction calculations for toric 
surfaces 

Subroutines Progm 

led into the Data General Eclipse 

CLOSE 

LENR, LENW 

Subroutine for closing a file 

FOPEN 

LENR, LENW 

Subroutine for opening a file 


9 



Table 5. Walk-Through of ICARUS Progran. 


To log onto ths tyBUm: 

USER NAME: 

SLfJ [saaple user inputs are indicated in this special type face] 
PASSWORD: 

STEVE [this can be changed] 

Aftor muth printing, ubbt may axeauta the progran by typing: 
t1)0 ICARUS 15 [this is systra dependent] 

IS THIS A NEW LENS? [Y or N] 

Y 

WHAT IS ITS NAME? [up to 7 alphanuaeric characters] 

LENSID [for exanple] 

If thie were not a new lens, the identifier etatment would oauee 
etdbroutine LEIfF (tene read) to be oalled; thie routine reads in 
stored lene data f^om memory. 

The following oomee from eubroutine STORE: 

'' ENTER OBJECT TYPE [the choices are: SUN (specifies sun as the object), 
PNT (point source), LIN (specifies an object defined in texns of 
its linear dimensions, ANG (specifies an object defined in terms of 
angular subtense) ; entering LIN or ANG triggers the following 
interactions] 


OBJECT DISTANCE IS: 
lOD [or whatever] 

LINEAR FIELD(S) « 

10 5 [height in y and x directions, respectively] 


OBJEa DISTANCE IS: 
lODD [ or whatever] 

ANGULAR FIELD(S) - 

5 □ [subtense in y and x directions, respectively; angles are 
in degrees] 


“I 


SUN 
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1 st surface 


ENTER ENTRANCE PUPIL DIAMETER [ctnrrmtly only circular q>tica can ba 
handled; the mtrance pupil diaaeter for aost aolar collectors will 
be the collector diaaeter] 

10 

ENTER DIUCNSIONS [M4 or IN or M or CM or FT, etc.; this could be changed 
to ENTER UNITS. Currently, this is of no use except perhaps to help 
the designer rearaber; in the future, prograa additims could aake 
use of this feature to scale output] 

DO YOU WISH TO SPECIFY RADII (R) OR CURVATURES (C)7 [R or C: since the 
prograa works in curvatures, if R is entered, all curve inputs will 
be divided into 1 to obtain curvatures; this pertains to surface data, 
toric curvatures, and Fresnel curvature] 

HON MANY LENS SURFACES ARE THERE? [generally, the answer will be 1, but 
the program can handle several; this ntjaber does not include object 
or ia^ge surfaces; the image is assumed to be plane and normal to the 
optical axis] 

S 

ENTER SURFACE TYPE [options are as follows] 

C(M: a conic section (ellipse, sphere, parabola, hyperbola) 

FNL: a normal Fresnel lens; aspheres are permitted 

CFN: curved Fresnel, a Fresnel lens that does not lie on a plane; 

only spherical base curves are permitted 
TOR: a toroid, a lens with different curvatures in the x and y 
directions 

TFN; toric Fresnel, a Fresnel with differing x and y curvatures 
CTF: curved toric Fresnel, a toric Fresnel that does not lie on a 
plane; (mly spherical base curves are permitted 

CON 

ENTER SURFACE DATA [respectively R, CC, TH, N^, N 2 , and N 3 , defined as 
follows; as in input elsewhere, numbers are separated by a space] 

R: radius of curvature or curvature of surface 

CC: conic constant: CC*0 sphere, CC*-1 paraboloid, CC > 0 -*• 

oblate spheroid, CC<-1 -»■ hyperboloid, 0>CC>-1 -► ellipsoid 

IH distance to vertex of next surface 

Ni,N 2 ,N 3 indices of refraction of three wavelengths 

SQQ >5 3.0 US3 ItSLA 1.5?7 

ENTER OECENTERS AND TILTS [respectively XDEC, YDEC (x- and y-decenters 
in coBBon units), a, 0 , y (tilts in degrees): these are currently not 
operational, Init so<m to follow are prograa additions that will in- 
corporate them; since they are currently disregarded by the program, 

0 is a good entry] 

□ soon 


ENTER RECENTERS AND BACmLTS [respectively RECx* REC2, a', Y*; 

again » these are not currently operational] 

Q □ □ □ 0 

ENTER SURFACE TYPE 
OF 

ENTER SURFACE DATA 

□ □ IQD 1 1 

[first two entries negated by choice of surface type; if a nontoric 
is requested in preceding line, then the R appears as the first 
entry here] 

ENTER DECENTERS AND TILTS 

□ □ □ □ □ 

ENTER RECENTERS AND BACKTILTS [same sign as original decenter or tilt 
cancels the former] 

□ 5 □ □ □ 

ENTER TORIC CURVATURES [CVX, CVY (x curve and y curve); this question 
appears only if a toroid has been specified] 

mi Plt TI [curvatures (C) or radii (R) , as previously specified] 

ENTER CURVATURE OF FRESNEL [base curve; appears if Fresnel is specified] 

“IbD [curvature (C) or radius (R) , as previously specified] 

Subroutine LENW (lens write) is now called; it creates a file called by 
the name the user specifies, and writes the lens to it. Should the r^ed 
arise, -this file may he accessed for modification without using the program. 

Return to main program: 

WOULD YOU LIKE TO REVIEW YOUR LENS? [Y or N; Y (yes) triggers subrou- 
tine LENT, which displays the tabulated data on the lens; this step 
is used to see whether any errors have occurred or whether the de- 
signer wants to change anything] 

Y 

Subroutine LEST (lens type): displays the current lens on the termi- 
nal in tabular fom. 

WOULD YOU LIKE TO MAKE ANY CHANGES? [Y or N; Y triggers subroutine 
CHANGE, which will then ask for the changes you'd like to make] 

Y 

Subroutine CEJ^GE : interacts with the user to make individual (flanges 
in the lens without having to retype the entire lens; all user inputs 
are as in subroutine STORE. 
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NHAT WOULD YOU LIICE TO DO NOW? [options are RAY, LENT, CHANGE, PARAX, 
RED, NEW, (^IT; handled by subroutines described below] 


I RAY 


1 


Subroutine RAY: traaee a specified ray through the lene. RAY in turn 
uses tuo subroutines, TORT and SYRM, to do the refraction catauZa- 
Hone for torio surfaces and axtatiy syimetrioal surfaces, respec- 
tively. Interaction is as follaoe: 

WHIOi COLOR IS TOE RAY, 1, 2, or 3? [1, 2, or 3: refers to Nj. N 2 , 

and Nj] 

ENTER FRACTIONAL OBJECT HEIGHTS [FOBY, FOBX: the nonualized (full 
field ■ 1) height on the object from which the ray is to start; ex- 
cluded if the object is SUN] 

ENTER FRACTIONAL PUPIL COORDINATES [RHOY, RHOX, the normalized (half 
of the entrance pupil diameter * 1) ray height on the entrance pupil; 
can be any coordinates, not the 0.7, 1.0 rays, which are built in 
when RED is required. The output is as follows: 




RAY COORDINATES AT IMAGE PLANE: 
R = (X2 + y2)^ 


X = 


Y = 


R = 


, 

I l-ENT I 

I Subroutine LENT (lens type): displays the current lens on the tevni- I 
' nal in tabular form. ' 


CHANGE 1 

I I 

I Subroutine CHANGE: interacts with the user to mike individual changes 1 

I in the lens without having to retype the entire lene; described < 

I above. 

I I 

I 1 

I PARAX I 

I Subroutine PARAX : does a paraxial ray trace that enables the user to 
I locate the Gaussian image plane, and in general to do a rough lens 
I design. One question is asked: 

• IS THIS IN COLOR 1, 2, or 3? [1, 2, or 3: determines whether Nj, N 2 , | 

* or Nj is used for each surface in the trace. Output is in the I 

I following form: I 

I SURF AXIAL Y CHIEF Y AXIAL U CHIEF U ] 

Q ***** , , . . 

I I ... . 
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I Subroutine RED: oalaulatee the radial energy distribution in the 
I image plane. The objeot is aUx^e assumed to be the sun. Inter- 
I action is as follcMs: 


• HOW MANY COLORS WOULD YOU LIKE? [1, 2, or 3: If 1 is specified, 

I index Nj is used and 96 rays are traced through the system. If 2 

I is specified, 96 additional rays are traced using N 2 at each sur- 

I face. If 3 is specified, 96 rays are traced using N 3 in addition 

I to those with Nj and N^. To do this, subroutine RAY is used. 

I The radial coordinate in the image plane is then stored, and sub- 

I routine SORT is called. This subroutine arranges the points in 

, order of increasing radius, and puts them into 10 bins. Each bin 
is now assumed to contain 10 % of the energy entering the system.] 

I Output is as folloQSy where = the radius containing successive 

I percentages of ray intercepts: 


PERCENT ENERGY 

RADIUS 

10 

ri 

20 


30 

T3 

100 

^10 


L I 


idiere SURF is the surface number (0 ■ <^ject; the last surface is 
the image surface) , AXIAL Y represents the height of the ray that 
passes through the center of the object and the edge of the 
entrance laipil, CHIEF Y is the height of the ray that passes 
throu^ the edge of the object and the center of the entrance 
pupil, AXIAL U is the angle the axial ray makes with the optical 
axis prior to refraction at specified surface, and CHIEF U is 
the angle of the chief ray; angles are measured counterclockwise 
from the axis] 


NCU [start on a new lens (or a different old one)] 


' llUIT [terminates the program; to restart it, type MXEQ ICARUS 15 
I (as before)] 

IHAT NEXT? [options are any of the above] 

auiT 

To log off Eclipse: 

BYE 


14 



original page is 

OF POOR QUALITY 


Table 4. User's Notes on Prograa ICARUS. 


1. The program takes 12 points to describe the sun, eight points on 
the lens in up to three wavelengths specified by three indices of 
refraction. This set then specifies 288 rays. If more points on the 
lens are desired, it will require recompiling the program and will in*- 
volve considerably more coiqjutation time (unless a single wavelength is 
desired, thus permitting 24 points on the lens). 

The pattern of points on the sun is a«^ shown below. They are dis- 
tributed in rings of equal area, each ring being weighted for a cosine 
solar limb-darkening effect. In reality, the power associated with limb 
darkening should be less than unity for closer approximations to the 
true sun, but for simplicity we have assumed a cos l.Ox relationship. 



Pattern of points on the sun distrilaited 
in weighted equal -area rings. 


2. The choice of surface type encompasses plane and plane tilted (as in 
the faceted Fresnel panel concept). 
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Note that the Fresnel surface has no finite grooves in this analyt- 
ical program. The proper tilt of the Fresnel element at the point where 
the entering solar rays reach the lens surface is automatically calcu- 
lated as discussed belw. 

The entrance pupil is in a plane normal to the vertex of the Fresnel 
base curve, as shcnm below. The present program takes eight rays: two 
radii (r » 0.7 and 1.0 the lens diameter at 6 « each of 0, 90, 180, and 
270 degrees. No central ray is traced. 


FROH SUN 



FRESNEL BASE 
CURVE (SPHERE) 


Definition of the pupil plane and Fresnel 
base (substrate) cuirve. 


The rays from the finite sun are directed to a single point (r,0). 
When the rays are transferred to the specified Fresnel base curve, the 
r,0 values will be slightly different because the rays from the disc of 
the sun are tilted with regard to the normal at the pupil surface. The 
same effect arises from either a tracking error or a deliberate tilt. 

The program calculates the normal at its intersection point of the 
Fresnel surface. The equivalent cur/ature of the Fresnel optical 
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surface is calculated, as illustrated below. The equivalent Fresnel ele- 
ment surfaces are always considered as curved elements, curved to yield 
a point focus at the specified TH value given for the final optical sur- 
face (if more than one is involved). 



BASE CUKVE 


Illustration of the difference between the 
base (substrate) curvature and the equivalent 
element curvature. 


The "proper tilt" of the Fresnel element has a subtle definition, 
as illustrated in the next figure. The equivalent Fresnel curvature has 
a constant value of RD (or CM) and CC, but is shifted by an amount A so 
that the slope at r brings the reflected ray at r to the focus. In 
other words, the program shifts the conic. 

Tbe radivis and curvature are defined by the focal length F or the 
power P, by 

CV « P/2 

P = 1/F 

RD equals 2F for a reflector, or 

1/F * (n' -n) (l/rj - l/r 2 ) ) 

I for a lens. 

P * (n* -n) (CVi - CV 2 ) ) 
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Geometry illustrating the shift in the 
effective conic surface defining facet angles. 


There is a strategy for finding the best value of CC for the equiv- 
alent facets, depending on the base surface used. For example, the best 
value of CC for a plane reflective Fresnel mirror is a hyperboloid, 

CC * -3.0. The strategy is to try a CC and examine the ray distribution 
in the focal surface, then change CC to find the best distribution. 


The curvature is the inverse of the radius of curvature, 
constant is: 


hyperboloid 

paraboloid 

allipse 

sphere 


CC < -1 
CC * -1 
-1 < CC < 0 
CC = 0 


The conic 


3. To find a better value of TH for the focus, take the output value 
of the y intercept with the TH plane and paraxial angle and calculate. 
Calculate an unproved TH using a pocket calculator. 



ATH » y* 
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Table S. Source Code Listing of Prograa ICARUS. 


PAIJUCBTBt «IAY-9(.ir-U,Or-Lf,1ir*2 
COMNOK/IUMI/ailMlI , XRZC ,LO , 1C .KOL 
CONMOR/OD/OBJ , DIM. KPD . LS .Junr .MX . unr . UK, TBf 
COMNON/SD/STmC If) ,CV( If) ,CC( If) .TH( If) ,IK 11,3) 

CONNOM/TD/DSC( If . 2 ) . TXLT( If . 3 ) . HDICl If . 2 ) . fTILT( If , 3 ) ,CYir( If) ,C7X( If) 
COMNOM/TD/CVF( If) 

COMM01)/ltaY8/rOBY( mtXY) , rOBX( MWiY) ,XTA( MAY) . YTA( IDtAY) . BA( HRAY) ,IIAfi( If) 
COMNOII/RAY8/R( NRAY, 3 ) 

REAL H.LMX.LMY 
LOGICAL LG 

If TYPE"XS THIS A DEW LENS?" 

READ(Xr.9ff) XAGE 
TYPE"WBAT X8 ITS NAME?" 

REAO(XF,9fS) ANEME 

XF( XAGE.EQ.*Y") CALL STORE 

XF( XAGE.EQ."N") CALL LENR 

TYPE"WOULD YOU LIKE TO REVIEW YOUR LENS?" 

READ(XF,9ff) LR 

IF( LR.EQ."Y") CALL LENT 

TYPE"WOULD YOU LIKE TO MAKE ANY CHANGES?" 

REAO(XF,9ff) XCQ 

1F(ICQ.EQ."Y") CALL CHANGE 

TYPE "WHAT WOULD YOU LIKE TO DO NOW?" 

R£AD( IF,9f5) WORK 

2f IF( WORK. EQ. "RAY") GOTO Iff 

IF{ WORK. EQ. "LENT") GOTO Ilf 
IF( WORK.EQ."PARAX") GOTO 12f 
IF( WORK. EQ. "RED") GOTO 13f 
XF( WORK. EQ. "CHANGE") GOTO 15f 
IF( WORK. EQ. "NEW") GOTO 3f 
XF( WORK. EQ. "QUIT") GOTO 4f 
TYPE"SORRY. I DON'T KNOW HOW." 

25 TYPE "WHAT NEXT?" 

READ(IF,9f5) WORK 
GOTO 2f 

Iff L6-.TRUE. 

CALL RAY 
GOTO 25 

Ilf CALL LENT 

GOTO 25 

12f CALL PARAX 

GOTO 25 

13f CALL RED 

GOTO 25 

15f CALL CHANGE 

GOTO 25 
3f GOTO 10 

4f CONTINUE 

9ff FORMAT( 81 ) 

905 FORMAT(87) 

END 

!A?lMETER*NRAY-f6,lF-ll,OP-lf,NF-2 

COMMON/ NAME /ANAME , IREC ,LG , 1C , NCOL 

COMMON /OD /OB J , D IM , EPD . LS .ANY , ANX , LNY . LNX , THf 

COMMON/ SD/8UR( 10) ,CV( 10) .CC( 10) ,TH( lf),N( 11.3) 

COMMOM/TD/DEC( 10.2). TXLT( 10.3). RDK( If. 2 ) .RTILTt If. 3 ) .CVY( If) .CVX( If) 
COM»«)N/TD/CVF< If) 

COMMON/ RAYS / FOBY( NRAY ) . FOBX( NRAY ) . XTA( NRAY ) . YTA( NRAY ) . RA( NRAY ) . RAD( If ) 
COM)«3N/RAY8/R( NRAY . 3 ) 

REAL N, LNY, LNX 
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Table 6. Source Code Listing of Prograa ICARI^ (Continued). 


TYPE"Birm OBJECT TYPE" 

EEE0(XF,999) OBJ 
1F(0BJ.BQ."SQM") THf-l.JE^lE 
IF(OEJ.B0.”8UH") GOTO If 
ACCEPT"OBJECT DXSTAECE IS ",THf 

XF(OBJ.EO."AEQ”} ACCEPT "AEGmAB FXELIMB) - ".AEY.AIIX 
Unr-THf «TAir( amy/57 . 296 ) 

XJOC-THf *TAM( AMX/57 . 296 ) 

XF(OBJ.BO."AMG") GOTO If 

XF(OBJ.R}."LXM") ACCEPT "LXMEAB FIELOCS) - ".XJnr.XJDC 
AMY-ATAM( LMY/TBf ) * 5 7 . 29 6 
AMX-ATAM( LMX/THf ) *57 . 296 

If ACCEPT "EMTER EMTRAMCE PUPXL DXMIETER "«BPD 

TYPE "EMTER DIMEM8X0MS" 

REAO(XF,999) DXM 

TYPE "DO YOU WX5H TO SPECIFY RADIX (R) OR CVRVATURE8 (C)?" 

READ<XF,998} IRQ 

ACCEPT "BOW MANY LEMS SURFACES ARE TBBRB7",L8 
DO 2f X-1«LS 

TYPE "EMTER SURFACE TYPE" 

READ(XF.999) SUR( X) 

ACCEPT "ENTER SURFACE DATAs ",CV( X) ,CC( X) ,TR( X) ,M( X,l) ,M( X,2) ,M( X,3) 
ACCEPT"EMTER DECENTERS AND TILTS " ,( DEC(I , J) , J-1 , 2 ) ,( TILT( I , J) , J-1 , 3) 
ACCEPT"EMTER RECEMTERS AND BACKTXLTS ",( RDBC( I, J) . J-1,2) .( RTXLT( X, J) , 

1 J-1,3) 

1F( SUR( X) .EQ.''TOR".OR.8UR( X) .EQ."TFM".OR.SUR( I) .EO."CTF") ACCEPT 
1 "ENTER TORIC CURVATURES " ,CVX( I ) ,CVY( I ) 

XF( 8UR( I) .EO."CFN".OR.8UR( I) .EO."CTF") ACCEPT"EMTER CURVATURE OF 
1 FRESNEL ",CVF(I) 

2f CONTINUE 

XF( XRO.FQ. "C") GOTO 4f 
DO 3f J-l.LS 
XF( CV( J) .EQ.f) GOTO 22 
CV( J)-1/CV( J) 

22 1F(CVY( J).EQ.f) GOTO 24 

CVY( J)-1/CVY( J) 

24 XF( CVX( J) .EO.f ) GOTO 26 

CVX( J)-1/CVX( J) 

26 XF( CVF( J) .EO.f ) GOTO 3f 

CVF( J)-1/CVF( J) 

3f CONTINUE 

4f CALL LENW 

998 FORMAT(SI) 

999 FORNAT<83> 

RETURN 

END 

SUBROUTINE LENR 

MIILMifiR HRAY-96 , XF-11 ,OF-lf ,MF-2 
COMMON / NAME /AMAFiE , X REC , LG , XC , MCOL 
COMMON/OD/OBJ , DIM , EPD , L8 . ANY , AMX , LMY , LMX , THf 
COMIK>N/SD/SUR( If) ,CV( If) ,CC( If) ,TH( If) ,M( 11,3) 

COMMON/TD/OEC( If ,2),TILT( If ,3) ,RDEC( If ,2) ,RTILT( lf,3) ,CVY( If) ,CVX( If) 
COMMON/TD/CVF( If) 

COMMON/ RAYS /FOBY( MRAY) ,FOBX( MRAY) ,XTA( MRAY) ,YTA( MRAY) ,RA(MRAY) ,RAD( If) 
COMIK)N/RAYS/R( MRAY, 3 ) 

REAL N,LNX,LMY 
CALL FOPEM( MF.AMANE) 

READ(NF,999) AMAME 

READ(NF,998) OBJ ,DIM,EPD,L8 ,AMY,AMX,LMY,LRX,TBf 
DO If I-1,L8 

READ<NF,997) SUR( I),CV( X) ,CC( X) ,TH( I) ,N( 1 , 1) ,N( X , 2 ) ,M( 1 , 3 ) 
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UJU)(Mr,9t6) (DEC(X,J).J-l,2).(TXI.T(Z,J),J-l,3),(B9aC(X,J).J-l,2}, 

1 (RTXLT<X.J),J-1.3) 

RSAD(1IP,995) CVY( X) ,CVX( X) ,CVT( X) 

If COfTXMnE 

995 rORMJlTi iX,3Flf.6) 

996 FORtiaTC lX,lfFlf.6) 

997 FOSMXTC 1X,83.6F12.6) 

998 FORM&TC 1X,83 , 2X,83 , 4Flf .6 , 3Blf .4 ) 

999 F08mT(lX,84) 

CALL CL08B(IIF.XBR) 

RETURN 

END 

8UBROUTINB LEWW 

^ARAHETER HMY-96 , XF-ll,OF-If ,irF-2 
CONMON/tlAME/ANANE . XREC , LQ , ZC , MCOL 
CONNON/OD/OBJ , DXM , EPD , L8 . ANY , AHX , Lmr , LRX , TSf 
COMNON/8D/8UR( If) ,CV( If) ,CC( If) ,TB( If) ,N( 11,3) 

CONNON/TD/l)BC( If ,2) ,TILT( If ,3) ,RDEC( If ,2) ,RTZLT( If ,3) ,C7Y( If) ,CVX( If) 
COMMOH/TD/CVF( If) 

CON>K>N/RAYS/FOBY( NRAY) ,FOBX( NRAY) ,XTA( KRAY) ,YTA( NRAY) ,RA( KRAY) ,RAI>( If) 
COMMON/RAYS/R( NRAY, 3 ) 

REAL N.LNY.LNX 
CALL FOPEN(NF,ARAHE) 

WRITE(MF,999) AKAME 

VfRlTE(NF,998) OBJ,DZM.EPD,LS,ANY,ANX.LNY,LIIX,THf 
DO 10 I-1,LS 

Vn)ITE(MF,997) 8UR( I) ,CV( X) ,CC( Z ) ,TH( X ) .N( 1,1) ,N( X,2) ,M( X,3) 
WRITE(MF,996) ( DEC(I , J ) , J-1 , 2 ) ,( TILT( I , J ) , J-1 , 3 ) ,( RDECC I , J) ,J-1,2), 

1 ( RIILT(I,J) ,J-1,3) 

WRITE(NF,995) CVY( I ) ,CVX( X ) ,CVF( I) 

If CONTINUE 

995 FORMAT( 1X.3F10 .6) 

996 FORMAK IX.lfFlf .6) 

997 FORMAT( 1X,S3,€F12.6) 

998 FORMAT( 1X,S3 , 2X.S3 ,4F10. 6 . 3Elf .4 ) 

999 FORMATC 1X.S4) 

CALL CLOSE(NF,ZER) 

RETURN 

END 

SUBROUTINE LENT 

NR^« 9 6 , IF- 1 1 , OF- If , NF- 2 
COMMON/NAME/ANAME , IREC . LG , IC , NCOL 
COMHON/OD/OBJ , DIN, EPD , L8 , ANY, ANX ,LNY , LNX , THf 
COMMON/8D/SUR( 10) ,CV( 10) ,CC( If) ,TH( 10) ,N( 11,3) 

COMW}N/TD/DEC( If ,2) ,TILT( If ,3) ,RDBC( If ,2) ,RTILT( If ,3) ,CVY( If) ,CVX( If) 
COMK)N/TD/CVF( If) 

C0NN0N/RAY8/F0BY( NRAY) ,FOBX( NRAY) ,XTA( NRAY) ,YTA( NRAY) ,RA( NRAY) ,RAD( If) 
C0MN0N/RAY8/R( NRAY, 3 ) 

REAL N,LNY,LNX 
DIMENSION RD( If) 

WRITE( OF.999) AKAME 
WRITE( OF, 899) 

WRITE( OF,998) OBJ, DIN 
WRITE( OF, 898) 

WRITE(OF,997) EPD, AMY, ANX, LNY, LUX, THf 
WRITE( OF, 897) 

DO If I-l.LS 

IF(CV( I) .EO.f) RD(I)-l.fE4>lf 
IF(CV( I) .NE.f) RD( I)-1/CV( I) 

If CONTINUE 

DO 2f I-1,LS 
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lfSXTI(OP,99<) CV( IMID( X),CC(1) ,Y« X).1K X.D r»( X.2),IRI.3} 

2f COHTXRDE 

WXXTK OF, 896) 

X-f 

DO 4f X-1,L8 

8D1I-DBC( X , 1 ) ••2*DBC( X , 2 ) **2*TX1T( X , 1 ) •*2*TXLT( X , 2 ) X , 3 ) ••2 
XF(8DM.IQ.f) GOTO 25 

mXTE(OF,995) ( DEC( X , J) , J-1 ,2 ) ,( TXLT( X , J) , J-1 , 3 ) 

25 8UM-XDEC( X,1)**2+EDEC{ X,2>**2+EIX1.T{ I,1)**2+KTXLT( X,2)»*2+*TXLT( X,3)*«2 

XP(8DM.B3.f) GOTO 3f 

WXITE(Or,994) ( BDEC( X,J),J-1,2) ,(RTXLT( X,J),J-1,3) 

GOTO 40 
30 X-K+1 

40 COMTIEUE 

Xr(X.E0.L8) WRXTE(Or,895) 

K>0 

WRITE( or, 894) 

DO 60 X-1,L8 

XF( 8UR( X).llE.''CrN".AllD.8UR( I) .EE. "CTF" .AllD.8tni( X) .EE.*'TOR".AIiD.8l7X( X) 

1 .HE."TFH") GOTO 50 

WRIIE(OF,993) CVY( X) ,CVX( X) ,CVF( I) 

GOTO 60 
50 K-K<^1 

60 CONTINUE 

Xr(K.EQ.LS) WRXTE(Or,89S) 

894 FORMAT(2X,"Y CURVE", 2X,"X CURVE" , 2X, "FRNL CURVE") 

895 rORIUT( IX , "NONE" ) 

896 PORMAT(2X,"y DBCENTER" , 2X . "X DECENTER" . 5X, "ALPSR" , 5X, "BETA" , 5X, "GANMA" ) 

897 rORHAT( IX , "CURVATURE ", 3X , "RADIUS ", 3X , "CONIC ", 2X , "THICKNESS " , 2X , "N1 " , 4X , 

1 "N2",4X,"N3") 

898 rORMAT( 6X,"BPD", 7 X." ANGULAR FIELDS" ,4X, "LINEAR FIELDS" , 5X, "OBJ. DISTANCE 

1 ") 

899 rORMAT( IX, "OBJECT UNITS") 

993 rORMATC 1X,3F8.6) 

994 rORHATC 1X,5F11.4) 

995 FORMAT( 1X,5F11.4) 

996 rORMAT( IX. 2F10.6 , F6 . 3 ,F10.4 , 3F7 .4 ) 

997 FORMAT( 1X.5F10.4.E10.4) 

998 FORMAT( 3X,S3,5X,S3) 

999 FORMAT( 1X,87) 

RETURN 

END 


10 


SUB 

paK 


ROUTINE CHANGE 

RLUeTER VUy-9« 


IP-11, OF-10,NF-2 


COMMON/NANE/ANAME , XREC , LG , IC ,NCOL 
COKNON/OD/OBJ , DIM , EPD . L8 , ANY , ANX , LNY ,LNX , TH0 
C0MM0N/8D/8UR( 10) ,CV( 10) ,CC( 10) ,TH( 10) ,N( 11,3) 

COMNON/TD/DBC( 10,2) ,TILT( 10,3) ,RDBC( 10,2),RTXLT( 10,3) ,CVY( 10) ,CVX( 10) 
COMMON/TD/CVF( 10) 

COMMON/ RAYS /FOBY( NRAY) ,FOBX( NRAY) ,XTA( NRAY) ,YTA( NRAY) ,RA( NRAY) ,RAD( 10) 
COMMON/ RAYS /R( NRAY, 3 ) 

REAL N.LNY.LNX 
INTEGER RO 

TYPE"MOULD YOU LIKE TO CHANGE THE 0BJK;T7" 

READ<XF,999) XTO 

IF( ITO.EO."N") GOTO 10 

TYPE"ENTER NEW OBJECT" 

READ(IF,998) OBJ 

TYPE "WOULD YOU LIKE TO CHANGE THE UNITS?" 
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MA0(Zr,999) ITO 
Xf(ZTO.SO."M") GOTO 2f 
TYPfSlITBIl mW tJKlTS" 

>EAD(Zr,998) DXM 

2t TYPE"WOm.D YOU LIKE TO CBJUME OBJECT DXSTBMCE2" 

EEAD(Xr,999> XTO 
IT( ITa.lO."H**) GOTO 3f 
ACCBPT'^HBW OBJECT DX8TMCE • ".TH0 
3f TYPE”VK)ULD YOC LIKE TO CRXM6E EPD?* 

READ(Xr.999) XTO 
Xr( XTO.EO."ll”> GOTO 4E 
ACCEPT"EBTEB MEW EPD",EPD 
4f TYPE”WOOLD YOU LXKE TO CHAMOE FXELD SXBE?” 

EEADiZF,999> XTO 
XF( XTO.BO."M*) GOTO 6E 
XF(OBJ.EO."LIM”) GOTO 5f 
ACCEPT"AMGULAR FIELDS - ".AMY.AMX 
LMY-TBS*TAM( AMY/S7 . 296 ) 

LMX-TB0*TAM( AMX/57 .296 ) 

GOTO 60 

50 ACCEPT"LXMBAR FIELDS - ",LMY,LMX 

AMY-57 .296*ATAM( LMY/TM0) 

AMX-57 . 296*ATAM( LMX/TH0) 

60 TYPB"DO YOU WISH TO SPECIFY RADIX ( R) OR CURVATURES (C)?” 

READ (IF, 999) RO 

TYPE"DOES THE LEMS STILL HAVE",LS," SURFACES?" 

READ ( IF, 999) ITQ 
IF( ITO.BO."Y") GOTO 70 
ACCEPT"EMTER THE MUMBER OF SURFACES", LS 
70 DO 90 Z-1,LS 

TYPB"DO YOU WISH TO CHAIRSE SURFACE ",I,"?" 

READ(IF,999) ITQ 
IF (ITO.EO."M") GOTO 90 
TYPB"EMTER SURFACE TYPE" 

READ(IF,998) 8UR( X) 

ACCEPT"BMTER SURFACE DATA< ",CV( I) ,CC( I) ,TH( X) ,M( 1,1) ,M( 1,2) ,M( X.3) 
ACCEPT"EMTER DECEMTER DATAi " ,( DBC( I , J> , J-1 , 2 ) ,( TXLT( X, J) , J*l. 3 > 
ACCEPT’EMTER RECEMTER DATAi " ,( RDECC I , J) , J-1 , 2 ) .( RTILT( X , J) , J-1 , 3 ) 

IF( SUR' I).BO."TOR".OR.8UR( I) .EO. "TFM" .0R.8UR( 1) .EO."CTF") 

1 ACCEPT"EMTBR TORIC CURVATURES" ,CVY( I) ,CVX( X) 

XF( SUR( Z).EO."CFM".OR.SUR( I) .BQ."CTF") ACCBPT"FMTER FRML CURVE" ,CV?( X) 
ZF(RQ.EO."C") GOTO 90 
XF(CV( I) .BO.0) GOTO 82 
CV( I)-1/CV( X) 


82 

XF( CVY( I) .BQ.0) 
CVY( I)-1/CVY( X) 

GOTO 

84 

84 

IF( CVX( X).BO.0) 

GOTO 

86 


CVX(I)- 1/CVX(1) 


86 

XF( CVF( X) .BO.0) 
CVF( I)-1/CVF( I) 

GOTO 

90 

90 

CONTINUE 



100 

CALL LEMW 



998 

F08NAT( S3) 



999 

FORNATC SI) 




RETURM 

BMD 

III^fiPirl{9^6,XF-ll,OF-lf,MF-2.XJI-L8*2 

CONNOR /MAKE /AMAME , X RK , LG , XC . MCOL 
CONMOH/OD/OBJ . DIM, %VD . LS , AMY , AMX , LMY, LMX , TH0 
CONNOR/SD/SUR( 10) ,CV( 10) ,CC( 10) ,TB( 10) ,M( 11,3) 
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COIQIOa/TD/DIC( If ,2) ,TXLT( If ,3) ,IS1C( If ,2) «KTILT( If .3) ,CVY( if) .C?X( if) 
CONNOM/TO/CVr( If) 

coMMOM/RaY8/roBY( maY) ,roix( nav) ,xra( mtaY) .yrac »iaY) ,ia( ntaY) .naiK if) 
coNNOii/aaYS/BC HaaY, 3 > 

DXNBMSXOM Y( 12 ) ,CY( 12 ) ,U( 12 ) ,CV( 12 ) ,»8X( 12 ) 

■BAT, M.LMY.LMX 

aCCEPT *XS TKXB XM COLOE 1.2,01 37", J 

Y(f)-f 

CY(f)— imr 

a(f)«EPD/{ 2*TRf) 

CU(f)— LMY/TBf 
PBI(f)-f 
Y( D-EPD/2 
CY( l)-f 

m l)-{0(fJ-Y( 1)»(II( 1,J)-1)*CV(1))/||{1,J) 

cm i)-ca(f)/ii( i,j) 

PHK 1)-(II( 1,J)-1)*CV( 1) 

DO If X-2,L8«1 

Y( I)-Y( I-1)^TH( X-D^m X-1) 

CY( X)-CY( I-1)^TH( I-l)«*Cm X-1) 

XF(I.EQ.L8<»1) GOTO If 

PHK I)-{ H( I,J)-H( I-1,J))*CV(X) 

m !)-(»( I-1,J)*0( 1-1)-Y( I)*PHI( X))/H( X,J) 

cm I)-( M( l-l,j)*cm 1-1)-CY( I)*PHI( X) )/H( X,J) 

M( I,J)-aBS( N( X.J) ) 

If CONTINUE 

PHI(L8'^l)-f 

m LS«l)-m L8) 
cm L8>l)«cm L5) 

WRXTE(OP,»99> 

DO 20 I-l, 1,8^2 
K-I-1 

WRZTE(Or,998) X,Y<K) ,CY(K) .mX) ,CmK),PBX(K} 

2f CONTINUE 

998 PORNAK lX,I4,Srif .€) 

999 PORMAK lX,"8URP",2X."aXiaL Y",SX, "CHIEF Y",3X, "AXIAL D",3X,"CHX*F U", 

1 5X, "POWER") 

CA7.L LENR 

RETURN 

END 

III^iilSltif-96 , IP-11 ,OP-lf ,NP-2 
CO)ON)N/NAME/ANANS. IREC.LG, ZC.NCOL 
C0)0«0N/0D/0BJ , DIM, EPD , L8 .ANY .ANX . LNY , LNX , TBf 
C0MM0N/8D/8URC If) ,CV( If) ,CC( If) ,TH( If) ,N( 11,3) 

COMION/ TD/DZC( If . 2 ) . TILT( 10,3), RDEC( If , 2 ) , RTILT( If . 3 ) , CTK If) , C7X( If ) 
COM)«DN/TD/CVP( If) 

C0M>K)N/RAY8/P0By( NRAY) ,POBX( NRAY) ,XTA( BRAY) . YTA( NRAY) ,IA( NJtAY) ,RAD( If) 
C0M)«)N/RAy8/R( NRAY, 3 ) 

CONNOR / TRAC E / YT , XT . A , L , M , ND 
COMI«>N/PUP/ RHOY( NRAY ),RHOX( NRAY) 

REAL N,L,N,ND,MU,LNY,LNX 

LOGICAL LG.LG2 

IP( .NOT. LG) NR-RRAY 

ZP< .NOT. LG) GOTO 30 

POBX< 1 ) -f 

POBY( 1 ) -0 

NR-1 

ACCEPT "WHICH COLOR 18 THE RAY, 1,2 OR 37", TC 
a%f,IC)-l 

ZP(0BJ.E0."8UN”) GOTO 10 
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AcciPT "iMTit riACTXoiua rajicT nxQns",roBY( i)«poBX(i) 

QOTO 2f 
if LMY-f 

UlX-f 

junr-t 

MX*! 

2f ACCBPt *BMTn PIACTXOIIAL P9PXI COORD XMATBS * .RHOYC 1) .XBOX( 1) 

If 00 2ff x-i,mi 

X^2-. FALSI. 

IX-LllX*rOBX< I) 

IY-UiY*rOBY( X) 

YT-RBOY( X)*IPD/2 
XT«RHOX( X)*IPD/2 

A-SQRT( THf*»2*( HY-YT)*»2*( BX-XT)»«2) 

L-(HX-XT)/A 
M-( HY-YD/A 
MD-TBf/A 
DO iff J-1,L8 

XF(ST7R( J).BO."TOR".OR.Btnt( J).»3.*CTr".OX.8UR( J).B0.*TF1I") TORTCJ.LOX) 

Xr(8DR( J>.BQ."C0II‘'.0R.8UR( J).BQ."PML*.0R.8tJR( J).BQ."CFM-) CALL 8YMN(J.L62> 
XP(L62) GOTO 2ff 
ISf COBTIRUE 

RXY-SQRTl XT**2^YT**2 ) 

XF(LG) TYPE "RAY COOROXKATES AT IMAGE PLAMIi X- ",XT." Y- ".YT," R- " 

1 .RXY 

IF(LO) GOTO 2ff 
XTA( I ) -XT 
YTA( I ) -YT 
RA( X>-RXY 
2ff CORTXMUE 

RETURN 
END 


CONNOM/MAME/AIIAMX , XREC , LG , XC , MCOL 
COMMOM/OD/OBJ .DIM, BPD.LS .AMY.AMX .LRY.LMX.Tllf 
C0MM0M/8D/8DR( if) ,CV( If) ,CC( If) ,TH( If) .R( 11.3) 

C010I0M/TD/DEC( If ,2) .TXLT( If ,3) ,RDBC( If .2) ,RTXLT( lf.3) .CYY( If) .C7X( If) 
CONMON/TD/CVF( If) 

CONNOR /RAYS /rOBY( RRAY) .POBX( RRAY) .XTA( RRAY) .YTA( RRAY) .RA( RRAY) .RAD( If) 
COKMN/RAYS/RC RRAY . 3 ) 

COMNOR/TRACE/YT . XT , A . L .N, RD 
COM»K)R/PUP/ RHOY( RRAY) .RHOX( RRAY) 

REAL R.L.N.RD.NU.LRY.LRX 
LOGICAL LG.LG2 
Y-Y«DEC( J , 1 ) >RDEC( J , 1 ) 

X-XOBCC J . 2 ) -ROEC( J , 2 ) 

R-8XR(A8XR(R)«TXLT( J.D/59.2IS78) 

M-8XR{ A8XR< M) ^TXLT( J . 2 ) /S7 . 29578 ) 

L-8XR( A8XR( L) <»TXLT( J. 3 ) /57 . 29578>RTXLT( J, 3 ) /57 .29578) 

N-8XR( A8XR<N)-tTILT(J. 2) 757.29578) 

R-SXR( A8XR( R)-RTXLT(J, 1)757 29578) 

Xr (8UR( J).BC."CTF"> GOTO If 
BX-CVX< J)*XT"*2 
BV-CVY{ J)*YT**2 
CX-BD-CVX( J ) *L"XT 
CY-RD-CFY( J)*M"YT 
DSX-CX««2-(CVX( J)*XT)*«2 
OS Y-C Y* • 2 - { C VY( J ) • YT ) • • 2 
GOTO 2f 
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1# BX-CVP( J)*{XT**2*YT**2) 

BY-BX 

CX-HD-CVr( J)*( L*XI*M*YT) 

CY"CX 

DSX-CX**2-CVP{ J)*BX 
D8Y-DSX 

2t IF(ABS(DSX) .LE.IE-IB) OSX-B 

XF(ABS(DSY) .LE.lE-10) DSY-0 
DX-8QRT( D8X) 

DY-SQHT( D8Y) 

A-(BX/(CX+DX)*BY/(CY+DY) )/2 

IF(8UR( J) .BQ."TFN") A-0 

XS-XT+A*L 

YS-YT+A*M 

8S-A*ND 

CSIX-DX 

CSIY-DY 

mj-H( J-1,IC)/ABS(H( J,IC) ) 

DEXS-1-KU'»*2*( 1-CSIX**2) 

DEYS-1-K0**2*{ 1-CS1Y**2) 

XF(DEXS.LE.0.OR.DEYS.LE.0) GOTO 70 

DEX-SQRT( DEXS ) -HU*CSIX 

DEY-SQRT( DEYS ) -MU*CS I Y 

ALPHA— CVX(J)*XS 

BETA— CVY( J)*YS 

L-HD*L+DEX*ALPHA 

M-MU*M+DEY»BETA 

MD-SQRT( 1-L**2-M**2) 

XT-XS+L-( TH( J)-ZS)/ND 
YT-YS*M*{ TH( J ) -ZS ) /ND 

A-SQRT( ( TH( J)-ZS)’»*2*( YT-YS ) **2*( XT-X8 > **2 ) 

GOTO 95 

70 XF(LG) TYPE"TOTAL XNTERKAL REFLECTXON AT 8GRFACE’',J 

GOTO 90 

80 XF(L6) TYPE”RAY MI8SED 8URFACE",J 

90 LG2-.TRDE. 

95 RETURH 

END 

SUBROUTINE SYMM( J.LG2) 

PARAHEfEftSfiAV-96, IP-11, OP- 10, NT-2 
COMMON / NAME / ANAME , I REC , LG , 1C . NCOL 
COMMON/OD/OBJ , DIM, EPD , LS , ANY , ANX , LNY , LNX , TH0 
COMMON/SD/SUR( 10) ,CV( 10) ,CC( 10) ,TH( 10) ,N( 11,3) 

COMI«)B/TD/DEC( 10,2) ,TILT< 10,3) ,RDEC( 10,2) ,RTILT( 10,3) ,CVY( 10) ,CVX( 10) 
COMMON/TD/CVF( 10) 

COMMON/ RAYS /FOBY( NRAY) ,FOBX( NRAY) ,XTA( NRAY) ,YTA( NRAY) ,RA( NRAY) ,RAD( 10) 
COMI«)N/RAYS/R( NRAY, 3) 

COMMON/TRACE/ YT , XT , A , L ,M, ND 
COMMON/PUP/ RHOY( NRAY ) , RHOX( NRAY ) 

REAL N,L,M.ND,KU,LNY,LNX 
LOGICAL LG, LG 2 
Y-Y+DECC J , 1 ) -RDEC( J , 1 ) 

X-X^DEC( J , 2 ) -RDEC( J , 2 ) 

N-5IN( ASIN< N)*TILT( J,l)/57. 29578) 

M-SIN( ASIN< M)>TILT( J,2)/57.29578) 

L-SIN( ASIN< L)+TILT( J,3)/57.29578-RTXLT( J,3)/57.29578) 

M-SIN( AS1N(M)-RTILT( J,2)/57.29578) 

N-SIN( ASIN( N)-RTILT(J,l)/57 .29 578) 

XP( 8UR( J) .BO."CPN") GOTO 10 
B-CV( J)*(XT*"2*YT*-2) 
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Table 5. Source Code Listing of Prograa ICARUS (Continued). 


C-M)-CV( J ) •{ I,*XT+M*YT ) 

D8-C*C-CV( 1+CC{ J)*I©*IID) 

GOTO 20 

10 B-CVT( J)*(XT**2*YT**2) 

C-KD-CVrC J) •( L*XT+K*YT> 

DS-C*C-CVP( J)*B 

20 IF(ABS(DS).LE.1S-10) D8-0 

D-SQRT( OS) 

A-B/(C4-D) 

IP(OBJ.BQ."WIL") A-0 

XS-XT+A»L 

yS-YT+A*M 

8S-A*ND 

E2-l/{ 1-2*CV( J)«CC( J)*8S) 

E-SQRT( E2 ) 

CSI-D*E 

MD-N( J- 1 , IC ) /ABS( N( J , IC ) ) 

DELS-l-MO**2*( 1-CSI**2) 

IFCDELS.LT.0) GOTO 70 

DEL-SQRT( DELS) -NU*CSI 

ALPHA- -CV( J ) »E*XS 

BETA— CV( J)*E*YS 

6AMKA-( 1-CV( J)*(CC{ J)+1)*ZS)*E 

L-i«7»L+DEL»ALPHA 

M-MO*M+DEL*BETA 

MD-MU*MD ^DEL *6AKMA 

XT-XS+L»( TH( J ) -2S ) /ND 

YT-YS+M*( TH( J ) -ZS ) /ND 

A-SORT( { TH( J ) -ZS ) **2 ♦( YT-YS ) »*2 ♦( XT-XS )** 2 ) 

GOTO 95 

70 XF(LG) TyPE"TOTAL INTERNAL REFLBCTION AT SURFACE ",J 

GOTO 90 

80 IF(LG) TYPE "RAY MISSED SURFACE ".J 

90 L62-.TRUE. 

95 RETURN 

END 


SUBROUT IRE RED 

ra5SHETBinniK7-96 ,IF-11 ,OF-10,«fF-2 

COMMON/NAME/ANAME , IREC ,LG , IC ,MCOL 
COMMOH/OD/OBJ , DIM, EPD , LS , ANY , AMX . LNY , LNX . TH0 
COMMON/SD/SUR( 10) ,CV( 10) ,CC( 10) ,TH( 10) ,N( 11,3) 

COMMOR/TD/DEC( 10,2) ,TILT( 10,3) ,RDEC( 10, 2 ) ,RTILT( 10 , 3 ) ,CVY( 10) ,CVX( 10) 
COMMON / TD/C VP ( 10) 

COMMON/RAYS /FOBY( NRAY) ,FOBX( NRAY) ,XTA( NRAY) ,YTA( NRAY) ,RA( NRAY) ,RAD( 10) 
COMMON/ RAYS /R( NRAY, 3 ) 

COMMON/PUP/RHOY( NRAY) ,RHOX( NRAY) 

REAL N, LNY, LNX 
LOGICAL LG 

DIMENSION FY( 12) ,FX( 12) ,RY( 8) ,RX( 8) 

DATA FY/, 31,0, -.31,0, .57,0, -.57,0, . 89 ,0, - . 89 ,0/ 

DATA PX/0, .31,0, -.31,0, .57,0, -.57,0, .89,0, -.89/ 

DATA RY/. 7,0, -.7, 0,1,0, -1,0/ 

DATA RX/0, .7,0,-.7,0,l,0,-l/ 

L6-. FALSE. 

LHX-4.36E+7 
LNY-4.36E+7 
DO 20 1-1,12 
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Table 5. Source Code Listing of Progra* ICARUS (Qmtinued). 


DO IB J-1,8 
K-8*( l-D+J 
P0BY(K)-ry( I) 

FOBX(K)-PX( 1) 

RHOY(K)-RY( J) 

BHOX(K)-SX( J) 
lir CONTIiriTE 

20 CONTINUE 

aCCEPT"BOW MANY COLORS MOULD YOTT LIKE? *,HCOL 

DO 40 J-l.NCOL 

CALL RAY 

DO 30 I-I,HRAY 

R( X,J)-RA( I) 

30 CONTINUE 

40 CONTINUE 

CALL SORT 
WRITE( OF. 999) 

DO 50 1-1,10 
J-10*I 

WRITE(OF.998) J,RAD(I) 

50 CONTINUE 

998 PORMAT( 7X.I3.7X.F10.4) 

999 FORMATC IX. "PERCENT ENERGY”, 5X. "RADIUS") 

RETURN 

END 

SUBROUTINE SORT 

PARAMETER NRAY-96 , IF-11 ,OF-10,NF-2 
COMMON/NAME /ANAME , IREC . LG , IC , HCOL 
COMMON/OD/OBJ .DIM, EPD ,LS , ANY , AHX .LNY ,LNX , TR0 
COMMON/SD/SUR( 10) ,CV( 10) ,CC( 10) ,TH( 10) ,N( 11,3) 

COMMON/TD/DEC( 10,2) ,TILT( 10,3) ,RDEC( 10,2) ,RTILT( 10,3) ,CVY( 10) ,CVX( 10) 
COMMON /TD/CVP( 10) 

COMMON /RAYS /FOBYC NRAY) ,FOBX( NRAY) ,XTA( NRAY) ,YTA( NRAY) ,RA( HRAY) ,RAD( 10) 
COMW3N/RAYS/R( NRAY, 3 ) 

NRA-NRAY"NCOL 
DIMENSION A( 288) 

REAL N.LNY.LNX 
DO 60 I-1,MRA 
Al I ) -0 

DO 20 K-l.NCOL 
DO 10 J-l.NRAY 
A( I)-AMAX1< A( I) ,R( J,R)) 

10 CONTINUE 

20 CONTINUE 

DO 40 K-l.NCOL 

DO 30 J-1,MRAY 

IF( A( I) .NE.R( J,K)) GOTO 30 

R( J,K)-0 

GOTO 60 

30 CONTINUE 

40 CONTINUE 

60 CONTINUE 

M-1FIX( NRA/10-».5) 

DO 70 1-1,9 

RAD( I)-A( NRA-*-l-M"I) 

70 CONTINUE 

RAD( 10) -A( 1) 

RETURN 

END 
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APPENDIX B 

EXAMPLES OF LENSES ‘INROUGH WHICH RAYS CAN BE TRACED 
1. Front Surface Reflective Fresnel 



2. Rear Surface Reflective Fresnel 
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3. Rear Surface Transaissive Fresnel 



Note that, when ’’Fresnel" is specified, the limitation of only a spheri- 
cal base curve enters. This limitation can be removed, but isn’t in 
this edition. The Fresnel elements, as specified by the "surface data,’ 
can be conic, toric, etc. 


4. Flat Transmissive Fresnel 


I 



Note w.iat a flat surface can be designated as a conic even though it is 
flat, and it will be treated as an asphere with the departure from the 
plane being the same as the departure of the corresponding conic from 
the sphere whose curvature it approximates at its vertex. 
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